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Abstract 

We  modified  the  lattice  fluid  equation-of- state  by  the  introducing  Debye-Huckel  equation.  A  thermodynamic  model  taking  into  account  the 
specific  interaction  and  ionic  strength  between  the  polymer  and  the  solvent  is  proposed.  The  proposed  model  successfully  predicts  the  vapor/ 
liquid  equilibria  (VLE)  of  solvents  and  the  solid  polymer  electrolyte  (SPE).  A  generalized  lattice  fluid  model  is  modified  to  describe  the 
change  of  water  activity  in  solid  polymer  electrolyte  (SPE)/water  systems.  The  calculated  activity  curves  using  the  proposed  model  agree 
remarkably  well  with  the  experimental  data. 
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1.  Introduction 

A  significant  time  generally  elapses  before  any  new 
technological  development  is  fully  exploited.  The  fuel  cell, 
first  demonstrated  by  Sir  William  Grove  in  1839,  has  taken 
longer  than  most,  despite  the  promise  of  clean  and  efficient 
power  generation.  However,  there  is  a  need  to  develop  more 
people  with  knowledge  of  fuel  cell  technology.  Today’s 
engineers  must  deal  with  increasingly  complex  and  non¬ 
ideal  pure  substances  and  mixtures,  and  with  operating 
conditions  involving  extremes  of  temperature,  pressure, 
concentration,  ionic  strength,  etc.  For  many  of  these  sub¬ 
stances  reliable  experimental  measurements  at  such  condi¬ 
tions  are  time  consuming  and  costly.  Using  the  simulation 
method  can  solve  these  problems. 

Perfluoro sulfonic  acid  membranes  are  the  objects  of 
intense  interest  due  to  their  potential  use  in  fuel  cells.  Since 
du  Pont  showed  that  the  presence  of  water  in  these  mem¬ 
branes  greatly  enhances  proton  transport,  considerable 
research  has  been  directed  toward  the  development  of 
SPE  having  high  ion  conductivity  as  well  as  mechanical, 
thermal  and  chemical  stability  [1,2].  One  important  factor 
that  can  affect  the  performance  of  fuel  cell  is  the  hydration  of 
SPEs  [3,4].  Membranes  can  neither  be  too  dry,  in  which  case 
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the  proton  conductivity  decreases  nor  too  wet,  in  which  case 
electrode  flooding  may  result.  To  understand  the  above- 
mentioned  phenomena  of  SPEs  for  the  fuel  cell  application, 
a  theoretical  study  is  desirable.  Investigating  the  manufac¬ 
turing  process  of  a  fuel  cell  in  common  use,  there  is  a 
pressure  imposed  on  the  SPE  to  improve  the  degree  of 
uptake  between  electrolyte  and  fuel. 

During  the  last  half-century,  there  have  been  various 
theoretical  and  experimental  approaches  to  investigate  the 
structures  and  thermodynamic  properties  of  polymers  and 
their  solvents.  The  first  description  of  statistical  thermody¬ 
namics  of  polymer  solutions  was  proposed  by  Flory  [5]  and 
Huggins  [6].  It  is  a  rigid  lattice  model  that  ignores  the  free 
volume  effect.  Since  then,  various  theories  have  been  devel¬ 
oped  that  allow  for  the  possibility  of  volume  changes. 
Prigigine  [7]  and  Eyring  et  al.  [8]  presented  a  cell  model 
theory  that  introduced  a  free  volume  by  assuming  a  hard- 
sphere  repulsive  potential  for  the  segments  of  the  polymer 
chain  moving  in  a  square- well  potential.  Flory  and  Flory 
et  al.  [9,10]  reported  a  modified  version  of  this  theory.  The 
main  difference  was  the  replacement  of  the  generalized 
Lennard- Jones  potential  by  a  van  der  Waals  type  potential. 
On  the  other  hand,  another  quantitative  theory  taking  the 
compressibility  of  both  components  into  account  has  been 
developed  by  Sanchez  and  Lacombe  [11-13].  This  is  similar 
to  the  Flory-Huggins  theory,  but  a  free  volume  term  is 
introduced  via  vacant  lattice  sites.  Based  on  this  lattice  fluid 
model,  Panayitou  [14]  derived  a  modified  molecular  theory 
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of  r-mer  fluid  mixtures  which  takes  into  account  the  differ¬ 
ence  of  the  molecular  shape.  Membranes  containing  ions 
with  multiple  charges  have  larger  effects  on  the  activity  of 
ions  than  membranes  containing  only  singly  charged  ions. 
To  express  membrane  ionic  contributions  in  a  way  that  takes 
this  into  account,  Debye  and  Hiickel  were  able  to  show  that 
in  solutions  the  activity  coefficient  of  an  ion  species  varies 
with  charge  number. 

The  purpose  of  this  study  is  to  modify  the  lattice  fluid 
model  to  the  polymer  solvent  systems  by  introducing  the 
Debye-Hiickel  model  as  a  function  of  the  polymer  ionic 
contribution.  We  examine  a  generalization  lattice  fluid  model 
for  water  activities  in  fuel  cell  systems.  If  a  membrane  is 
contacted  with  an  aqueous  solution,  the  membrane  will  take 
up  water  and  solutes.  When  the  membrane  with  its  dissolved 
water  and  solutes  is  considered  as  an  aqueous  solution,  the 
fixed  charge  groups  of  the  membrane  are  ordinary  ions  evenly 
distributed  in  the  aqueous  solution. 


2.  Model  development 


2.2.  Chemical  potential 


The  free  energy  of  mixing  is  a  minimum  when 
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(, 


T 


pz  -t-P-f-7^1n(l  -p)  +  ^1  --JP  )  =  0 
where 


p  = 


and 


1 


kT 


v  = 


T  =  —  = 


=  <^^1  +  4>2v  *2, 

i 

T  kT 


~  _  P 

P  *  1 

p * 


(4) 


(5) 


Two  theoretical  aspects  are  taken  into  account;  the  GLF 
model  and  Debye-Hiickel  theory.  We  assume  that  the 
solvent  is  a  component  in  order  to  take  account  of  a  solvent 
effect  of  the  free  energy  of  mixing.  We  assume  that  the  salt 
is  component  in  order  to  take  account  of  a  salt  effect  of  the 
free  energy  of  mixing. 

In  this  study,  all  expressions  for  the  thermodynamic 
model  framework  are  expressed  in  terms  of  a  binary  mixture 
of  N i  molecules  of  size  r{  and  N2  molecules  of  size  r2-  The 
r-mer  is  a  molecule  that  occupies  r  sites  on  a  lattice  of 
coordination  number  z. 

2.1.  The  free  energy  of  mixing 

In  the  generalization  lattice-fluid  (LF)  model,  free  energy 
per  molecule  G/rN  as  fE,  is  given  by 

Pv* 

fE  =  — pe*  +  —3 - TS  (1) 
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where  s*  is  the  mixing  interaction  energy. 


2.3.  Debye-Hiickel  theory 


To  take  into  account  the  ionic  contribution  between  the 
segments  of  molecules,  Debye  and  Hiickel  [15]  showed  that 
in  dilute  solutions  the  activity  coefficient  yt  of  an  ion  species 
i  with  a  charge  number  of  n  is  given  by 

logy,  =  -AzJ1/2  (6) 


where  I  is  the  ionic  strength  and 
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where  rasolv  is  the  mass  of  solvent  in  volume  V,  e  is  the  unit 
charge  1.621  x  10“19C  and  sr  is  the  relative  permittivity. 

At  higher  ionic  strength,  the  following  empirical  exten¬ 
sion  is  often  useful: 
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where  B  presents  the  dependencies  of  the  mean  ionic  activity 
coefficients  of  components  on  ionic  strength  for  the 
(2)  extended  Debye-Hiickel  theory. 


fi2  is  the  free  energy  parameter  that  replaces  the  pure 
energetic  parameter  e12.  It  includes  the  specific  interaction 
parameter  Ss/k  and  the  entropy  effect  generated  by  the 
specific  interaction. /i 2  is  defined  as  [16] 
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Subscript  1  refers  to  1 -component  (solvent)  and  subscript  2 
refers  to  2-component  (polymer). 


2.4.  The  activity 

The  activity  of  a  mixture  is  defined  by 
Pi  =  p°RT\nai  (9) 

Thus,  the  activity  at  is  simply  a  means  for  expressing 
the  chemical  potential  of  a  species  on  a  mixture.  In  general, 
ion  exchange  membranes  uptake  water  component.  Ther¬ 
modynamics  gives  this  chemical  phenomenon  to  activity. 
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In  a  binary  mixture,  the  activity  of  solvent  1  in  polymer  2 
is 
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3.  Results  and  discussion 

We  employed  the  modified  generalization  lattice  fluid 
model  to  predict  water  activities  of  perfluorosulfonic  acid 
polymer  electrolyte/water  systems.  The  advantage  of  this 
model  is  its  simplicity. 

Fig.  1  shows  basic  structures  of  Nation,  Aciplex  and 
Flemion.  In  SPE,  the  HS03  group  added  is  ionically  bonded, 
and  so  the  end  of  the  side  chain  is  actually  an  S03“  ion.  The 
result  of  the  presence  of  these  S03-  and  H+  ions  is  that  there 
is  a  strong  mutual  attraction  between  the  +  and  —  ions  from 
each  molecule.  Now,  a  key  property  of  sulfonic  acid  is  that  it 
is  highly  hydrophilic — it  attracts  water.  In  SPE,  this  means 
we  are  creating  hydrophilic  regions  within  a  generally 
hydrophobic  substance.  The  hydrophilic  regions  around 


—  (CF2CF2)x— CFCF2  — 

ocf2  cfocf2cf2  —  so3h 

(a)  CF3 


—  (CF2CF2)x  —  cfcf2  — 

(OCF2CF)mO(CF2)n— S03H 

(b)  CF3 

Fig.  1.  Repeating  unit  structures  of  perfluorinated  proton  exchange 
membrane  with  sulfonic  acid  functional  groups,  (a)  Nafion  (x  =  5-11)  (b) 
Aciplex  and  Flemion  (m  =  0  or  1,  n  =  2-5,  x  =  1.5-14). 


Water  collects 
around  the  clusters 
of  hydrophilic 
sulphonate  side 
chains 


Fig.  2.  Showing  the  structure  of  Nafion  type  membrane  materials.  Long 
chain  molecules  containing  hydrated  regions  around  the  sulfonated  side 
chains. 


Fig.  3.  Water  activities  for  the  Nafion  (E.I.  duPont  de  Nemours  and  Company)/water  system  at  353.15  K  Nafionll7  (EW  =  1100).  The  solid  line  is 
calculated  from  this  work.  The  other  line  is  for  the  original  generalization  lattice  fluid  model.  Solid  squares  are  experimental  data  [18]. 
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Fig.  4.  Water  activities  for  the  Aciplex  (Asahi  Chemical  Industry  Co.)/water  system  at  353.15  K.  Aciplexl2  (EW  =  1080).  The  solid  line  is  calculated  from 
this  work.  The  other  line  is  the  generalization  lattice  fluid  model.  Solid  squares  are  experimental  data  [18]. 


the  cluster  of  sulfonated  side  chains  can  lead  to  the  uptake 
of  large  quantities  of  water.  Fig.  2  shows  the  micro-phase 
separated  morphology  [17].  Fig.  3  presents  calculated  water 
activities  of  Nafionll7  (du  Pont,  EW  =  1100)  with  experi¬ 
mental  data  reported  by  Hinatsu  et  al.  [18]  We  compared 
original  GLF  with  modified  GLF  (including  ionic  contri¬ 
bution).  The  each  parameter  values  (subscript  1:  solvent, 
subscript  2:  polymer)  of  Nafion  117/water  system  are 
&\2  =  — 554.48  K,  Ss  =  — 14.53  K  and  ionic  parameter 
B  =  3.992.  These  parameters  are  determined  by  fitting 
experimental  data.  The  proposed  model  agrees  well  with 


all  regions.  However,  the  GLF  equation  overestimates  in  the 
low  water  content  region.  Fig.  4  presents  the  water  activities 
for  Aciplex  12  (Asahai  Glass,  EW  =  980).  The  each  para¬ 
meter  values  are  &\2  =  —566. 39  K,  ds  =  — 19.99  K,  and 
ionic  parameter  B  =  0.274,  respectively.  Fig.  5  shows  the 
activities  for  Flemion4  (Asahai  Chemical,  EW  =  890).  The 
parameter  values  are  s*l2  =  —  577.58  K,  3s  =  —  25.16  K,  and 
ionic  parameter  B  —  — 0.5836. 

The  ultimate  goal  of  this  work  is  to  improve  fuel  cell 
properties  such  as  methanol  crossover.  As  can  be  seen  in 
these  figures,  the  water  activities  increase  as  the  equivalent 


Fig.  5.  Water  activities  for  the  Flemion  (Asahai  Glass)/water  system  at  353.15  K  Flemion4  (EW  =  890).  The  solid  line  is  calculated  from  this  work.  The 
other  line  is  for  the  generalization  lattice  fluid.  Solid  squares  are  experimental  data  [18]. 
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weight  of  the  polymer  decreases,  or  in  other  words,  as  the 
number  of  exchange  sites  increases.  The  ionic  strength  plays 
a  deceptive  role  in  this  system. 

In  this  study,  we  determine  interaction  energy  parameters 
and  ionic  parameters  between  perfluorosulfonic  acid  polymer 
electrolyte  and  the  water  component.  The  results  obtained  by 
the  proposed  model  are  expected  to  provide  the  appropriate 
operating  conditions  for  fuel  cells.  They  are  important 
because  cell  performance  can  be  adversely  affected  if  the 
SPE  is  allowed  to  dry  out  (resulting  in  lower  membrane 
conductivity)  or  if  excess  water  is  allowed  then  electrode 
flooding  occurs  (resulting  in  lower  electrode  performance). 
To  characterize  most  common  polymer  electrolyte/fuel  sys¬ 
tems,  more  experimental  data  with  different  pressure  systems 
are  required. 

In  the  context  of  ion  exchange  membrane  fuel  cells,  the 
membrane  activity  is  not  the  only  issue.  From  the  scientific 
viewpoint  one  must  also  consider  aspects  such  as  water 
transport  through  the  membrane,  chemical  and  mechanical 
stability,  the  ability  to  support  electrochemistry,  and  the 
formation  of  low  resistance  interfaces  with  gas  diffusion 
electrodes.  From  a  technological  viewpoint,  processability, 
reproducibility,  and  cost  must  also  be  considered.  In  fuel 
cells,  membranes  operate  under  significant  mechanical  pres¬ 
sure  that  acts  in  opposition  to  the  osmotic  pressure.  Never¬ 
theless,  these  results  should  serve  to  help  in  understanding 
proton  conductivity  in  PEMs  and  assist  the  design  of  PEMs 
with  even  lower  ohmic  resistance. 

4.  Conclusion 

This  work  showed  that  activity  in  ion  exchange  mem¬ 
branes  is  determined  by  several  factors.  To  describe  the 
activity  these  of  systems,  we  modified  a  generalization 
lattice  fluid  model.  The  model  has  a  simplified  and  improved 
expression  for  the  free  energy  of  mixing  for  polymer/solvent 
systems  that  includes  the  combinatorial  and  vacant  entropy 
contribution,  the  van  der  Waals  energy  contribution,  and  the 


ionic  contribution.  We  compared  calculated  water  activities 
of  several  commercial  perfluorosulfonic  membranes  with 
experimental  sorption  data.  The  activity  came  under  the 
influence  of  the  ionic  contribution.  This  model  gives  good 
agreement  with  experimental  data. 
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